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Triangular symmetry stabilizes a novel non-Fermi-liquid phase in the three-impurity Kondo model 
with frustrating antiferromagnetic interactions between half-integer impurity spins. The phase arises 
without fine-tuning of couplings, and is stable against magnetic fields and particle-hole symmetry 
breaking. We find a conformal field theory describing this phase, verify it using the numerical renor- 
malization group, and extract various exact, universal low-energy properties. Signatures predicted in 
electrical transport may be testable in scanning tunneling microscopy or quantum-dot experiments. 
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The same many-body physics that is responsible for 
the Kondo screening of magnetic impurities in bulk met- 
als 0] produces resonances in tunneling through a quan- 
tum dot or an adatom on a metallic surface 
Greater experimental control over the latter settings al- 
lows systematic study of multiple- "impurity" configura- 
tions in which the Kondo effect competes with ordering 
of the local degrees of freedom ■ One goal that 
remains elusive is experimental realization of non- Fermi- 
liquid (NFL) behavior similar to that exhibited, ej^jby 
the two-impurity and two-channel Kondo models |7], |8| • 

A cluster of three antiferromagnetically coupled spins 
is of fundamental importance as the simplest example of 
frustration, a feature of many magnetic systems. Scan- 
ning tunneling microscopy (STM) has found two distinct 
types of compact Cr trimer on a gold surface || : "type-2" 
trimers show a sharp resonance of width (Kondo temper- 
ature) Tk ~ 50 K, whereas any Kondo effect for isolated 
Cr atoms and "type-1" trimers seems to have TV <C 7K. 
Attempts to explain this result via variational 9] , quan- 
tum Monte Carlo [10| , and perturbative renormalization- 
group (RG) ^3 treatments of a three-impurity Kondo 
model have reached opposing conclusions concerning the 
tri ang ular geometry of the type-2 trimers: equilateral 

[llj or isosceles Qj| • Interest is also developing in the 
interplay between Kondo physics and interdot quantum 
entanglement in triangular quantum-dot devices [l^. 

This Letter reports exact results on a frustrated phase 
of the three-impurity Kondo model exhibiting NFL be- 
havior distinctly different from that found in previ- 
ously studied models. We present a conformal field 
theory (CFT), deduced by comparison with numerical 
renormalization-group (NRG) results, showing that the 
NFL fixed point is stable against particle-hole-symmetry 
breaking (unlike its two- impurity counterpart), exchange 
anisotropy, and even an applied magnetic field (which 
destroys two-channel Kondo NFL behavior). This en- 
hanced stability compared to other NFL fixed points 
makes the frustrated phase an excellent candidate for ex- 



perimental realizations; indeed, it has been argued on the 
basis of weak-coupling RG to describe the low-energy 
physics of the type-2 Cr trimers in [j| . We make predic- 
tions for the conductance expected in STM experiments 
on trimers and in certain quantum-dot devices. 

Model. — We start with a Hamiltonian /Tband + H mt de- 
scribing a noninteracting conduction band coupled via 
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to spin-5 impurities Sj (J = 1,2,3) at the vertices rj 
of an equilateral triangle; ip a ij) annihilates an electron 
with spin a = ±| at r . We assume that the permu- 
tation group S3 that maps the set {r}} onto itself is a 
subgroup of the lattice symmetry group (as is the case, 
e.g., in The impurities couple to just six orthonor- 
mal combinations of conduction states, annihilated by 
operators iph,a e" ,2 ^' 1, ' 3 ^ Q (r :J -), where h = 0, ±1 

is the "helicity": under a 27r/3-rotation about the center 
of symmetry, a helicity-/i state is multiplied by e l27Th ^ 3 . 
The combined states of the three impurities can also be 
constructed to have well-defined helicities, in which case 
the Hamiltonian conserves total helicity (modulo 3) and 
is invariant under interchange of all helicity labels 1 and 
— 1. Then, Eq. Q) can be rewritten (l3| 

-ffint = [Too sqo + Jn(sii + sn)]-5o + [Joi(sbi + si ) 
+ Tiisii]-5i + [T i(sio + s i) + Tnsn] -Si, (2) 

where S h = £\ e^^Sj, s w = £ Q/3 ^ h ' a \B a ^ w ^ 
and 1 = — 1; Jhh' equals T times a non-negative fac- 
tor that depends on the impurity separation and the 
conduction-band dispersion, as well as h and h! [T^J. 

For 8 = 2 , the NRG shows 01 that over a large region 
of the parameter space of Eq. J2J) , the low-energy physics 
is governed by a "frustrated" fixed point at which the 
impurities are locked into the subspace of two doublets 
of combined spin Si mp = \, one each of helicity h = ±1. 
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Three spins of arbitrary half-integer S", coupled by an 
additional Hamiltonian term K J2i<j &i ' with K J , 
also lock into an Simp = \, h — ±1 subspace. Weak- 
coupling RG analysis of this augmented model, which 
for S — | provides a description of equilateral Cr trimers, 
is consistent with flow to the same fixed point; for S = 
|, moreover, the characteristic temperature Tk for this 
flow is found to greatly exceed the single-impurity Kondo 
scale, in agreement with the Cr-trimer experiments 0. 

In the frustrated phase, Jn in Eq. © scales to zero, 
Joo and Jn can be neglected, and particle-hole asymme- 
try is marginal Thus, we analyze the fixed point in 
a restricted Si mp = \ space, replacing Eq. J2J by 

H int = -V^M^^aT+^Tr^ + H.c] ■ S imp . (3) 

Here, and T z act on the electron helicity in the spin- 
1 representation of an "orbital-spin" SU^(2) ^5}, with 
matrix elements (T z ) h , h = h, (T+) 1>0 = (T+) _i = V2. 
The Pauli matrices a, 2S lmp = 2Sq and f\ mp act, respec- 
tively, on the electron spin, impurity spin, and impurity 
helicity, with 0f mp ) Ml = -h for h = ±1, (t+ 1p )_ m = 1. 

It is important to note that setting Jjj = enlarges the 
S 3 symmetry of Eq. J2J to a {/(') (1) symmetry in Eq. 
replacing total helicity (conserved only modulo 3) by a 
conserved quantity t z : the eigenvalue of i/^T 2 ^ + ^T z mp . 
Now, commutes with SU (s) (2) spin, orbital 
spin, and also with SU^{2) isospin defined by I z = 

CFT description. — We obtain exact analytical results 
for the frustrated fixed point using the boundary CFT 
approach to quantum impurity problems 16]. The op- 
erator iph,a is considered to act at the boundary x = 
of a one-dimensional space < x < I 0. The key is 
to find a "conformal embedding " (a decomposition of 
the bulk fermions iph,a(x) into products of "constituent" 
fields admitting a "fusion procedure" (generating a 
new conformally invariant boundary condition) that re- 
produces the fixed-point finite-size spectrum (FSS). We 
deduce this FSS by extending to higher accuracy the 
NRG results of 0] . From the fusion procedure, all uni- 
versal low-energy properties in the physical limit I — > oo 
can in principle be computed exactly. 

We first construct a conformal embedding of the free 
Dirac fermions iph,a{x) in which the hclicitics transform 
in the spin-1 representation of an SU^(2) "pseudospin" 
P, where P z = T z and P + has matrix elements in the 
helicity basis (P+)i, = -(T+) li0 , (P + )o,-i = (T+) ,-i. 
Unlike T defined above, P commutes with isospin I, The 
free-fermion FSS can be decomposed into products of 
SU( S \2) 3 x SZ7«(2) 3 x SU&> (2) 8 conformal towers [3 
as exemplified in Table [I] for boundary conditions that 
yield a nondegenerate ground state. Here, SU(2)k is a 
level-fc Kac-Moody CFT; see and references therein. 

Since Eq. © lowers the SU ( - P \2) symmetry of P ba nd 
to t/W(l) = C/ (t) (l), we analyze the frustrated fixed 



TABLE I: Finite-size spectrum of free fermions decomposed 
into products of spin, isospin, and pseudospin conformal tow- 
ers, labeled by s, i, and p, respectively. The subscript A gives 
each tower's contribution to the excitation energy jlfll ] . 



(s)a 


(*)a 


(p)a 


(0)o 


(0)0 


(0)o + (4)2 


(1)3/4 


(D3/4 


(0),, + (4) 2 


(2)3/20 


(5)3/20 


(l)l/5+ (3) 6 /5 


(1)2/5 


(1)2/5 


(l)l/5 + (3) 6 /5 


(0)o 


(1)2/5 


(2)3/5 


(1)2/5 


(0)o 


(2)3/5 


(2)3/20 


(3)3/4 


(2)3/5 


(f)3/4 


(2)3/20 


(2)3/5 



point using an embedding obtained from that above by 
decomposing 5C7 (p) (2]s D ;7 (t) (l) 8 x Z 6 , where Zg is a 
parafermionic CFT [fj. In any SU(2) k CFT, each pri- 
mary operator 0w , transforming in the spin-j represen- 
tation (j = 0, i, 1, fc/2), factors into a sum of products 
of a Zfr primary ip-^ and a U(l)k boson exponential |2d| : 

= ^ ^e^^P. (4) 

Setting k — 8 and j = p, we rewrite the spectrum in 
Tabledas products of SU^(2) 3 x SU®(2) 3 x E/W(l) 8 x 
Zg, conformal towers. These products also provide the 
operator spectrum for the noninteracting model. 

Using this conformal embedding, we are able to obtain 
the frustrated fixed-point FSS from the free-fermion FSS 
by applying a three-step fusion procedure: (I) fusion with 
the s = I primary operator in SU^(2) 3 , then (2) fusion 
with the V = \ primary operator in SU^ (2)g, then (3) 
fusion with the 4*^=2 primary operator in Zg |2*H . 

As an illustration, Table ITT1 lists all CFT states of en- 
ergy 0| E < 1 in the frustrated FSS for boundary con- 
ditions that yield a degenerate free-fermion ground state 
|22j , along with energies of NRG levels having the same 
(s, i, t z ) quantum numbers. To fix the overall NRG en- 
ergy scale, distorted by band discretization, we match the 
lowest excitation of the free-fermion spectrum to its CFT 
counterpart. Apart from small energy shifts (residual 
discretization effects), the CFT and NRG spectra agree 
perfectly. In fact, all 1810 CFT states with E < 1.8 have 
been compared with and match NRG levels . 

Applying the fusion procedure twice to the free- 
fermion FSS gives the complete and exact spectrum of 
boundary operators that can be added to the fixed-point 
Hamiltonian ^(|. This spectrum (see Table ITTTll exhibits 
"fractionalization" of charge, spin, and orbital degrees of 
freedom, as is typical of an NFL fixed point. Remarkably, 
it also exhibits full SU^{2) symmetry. The SU^{2) 8 
operator multiplets (last column of Table 1111(1 should 
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TABLE II: Finite-size spectrum at the frustrated fixed point. 
U^'(l) and Zg conformal towers are labeled by t z and (p,m), 
respectively. Each row represents a pair of states related by 
a change in the signs of t z and m. E is the CFT excitation 
energy and Snrg is the NRG energy computed for a band 
discretization parameter A = 3. See also Table [I] 



(s)a (i)a 



(p,m)A 



E E 



NRG 



(0)o 



(0)o 



(2)9/32 



(2 



(2)3/20 (1)3/20 (1)1/3 



2 



(0) o 

(1) 2/5 



(1)2/5 
(0)o 



(2)1/32 
(2)1/32 



(- 

\ 2 



(§) 



2)3/20 



•.2)3/20 



(1)9/32 



(2 



(0)o 
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(0)o 
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1/32 
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(2 
(2 
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(2)3/20 
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(2 
(I 

(h 
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9/32 
1/32 
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V 2 


+ |)l9/32 

-Dll/32 
3\ 
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0.6 
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(h 
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3 \ 
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(^ 
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0.7 


0.7043 


(1)2/5 


(0)o 


(2 
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(2 
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1/32 
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(I 
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9/32 
1/32 


(£ 
(2 


— 5)7/160 
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(I) 



25/32 



(3 



-3WI6O 



3\ 
-5j3/32 



0.1 0.1001 



3\ 
"2i3/32 

3\ 
-5J3/32 



0.2 0.2000 
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-5J7/160 



0.3 0.2996 



3\ 

-5J127/I6O 
'1)7/160 



0.5 0.4968 
0.5 0.5020 



+1) 



7/160 



0.8 0.8116 



again be decomposed using Eq. into U^(l)s x Z%. 

Boundary operators entering the effective low-energy 
Hamiltonian for the frustrated fixed point must respect 
the SU<i)(2) x SU^(2) x f/W(l) symmetry of the full 
Hamiltonian Such operators appear in the first row 
of Table ITTT1 Only (s,i,t z ,Z a ) = (0,0,0,(^0)1/5) is rel- 



TABLE III: Operator spectrum at the frustrated fixed point. 
A gives each factor's contribution to the scaling dimension. 
"2x" indicates two operators with the same p and A. 
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(0)o 


+ 2x[(l) 1 /5 + (2) 3 /5 + (3) 6 / 5 ] + (4)2 


(1)2/5 


(1)2/5 


(0)o 


+■ 2x[(l) 1/5 + (2) 3 /5 + (3) 6/5 ] + (4) 2 


(0)o 


(1)2/5 




(l)i/5 +2 x (2)3/5 + (3) e/5 
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evant (in the RG sense). It cannot appear because it 
is odd under the Z2 subgroup of S3: iph,a ~> ~i } -h,a, 
T^r, T^t„, which is representable as a 7r-rotation about 

liiljj imp ' ■*- 

the cc-axis in orbital-spin space |23j. The least-irrelevant 
operator also respecting this discrete Z2 symmetry of 
Eq. is the corresponding S{/( p )(2) 8 -descendant of 
dimension A = 1 + 1/5, which yields a correction-to- 
scaling exponent 1/5 in excellent agreement with the 
value 0.200 ± 0.002 observed in the NRG spectrum. 

Physical results. — We now present exact properties 
that can be deduced from the CFT description. Details, 
including analysis of the conditions required for observa- 
tion of these properties, will appear elsewhere [2^. 

(a) Fixed-point properties. — The frustrated fixed point 
has an irrational "ground-state degeneracy" [24| g = 
[|(5 + Vb)] 1 / 2 . Moreover, in a quantum-dot device of 
triangular symmetry where biases Vj in leads j = 1, 2,3 
produce in lead i a current Ii = 53 . Vj , the T = 
zero-bias conductance is Ga = 4e 2 /3/i. By contrast, the 
"isospin two-channel" regime [l3T |. in which Jjj domi- 
nates Eq. J2J, is unstable against particle- hole asymme- 
try and at low energy exhibits the Fermi-liquid behavior 
of the SU (4) fixed point of 0(a)], with g = 1 and (in the 
limit of small particle-hole asymmetry) Gu = 8e 2 /9h < 
4e 2 /3/i. The other stable fixed point of [13J, at which 
inter-impurity correlations are irrelevant and the stan- 
dard Kondo effect is recovered, has 3 = 1 and Ga = 0. 

(b) Differential conductance. — The leading irrelevant 
operator of dimension A = 1 + 1/5 governs many prop- 
erties near the fixed point. In particular, the differential 
tunneling conductance into the impurities from a metal- 
lic lead (e.g., an STM tip located symmetrically with re- 
spect to the impurities) in the regime k B T, \eV\ <C k B Tx 
(V being the bias voltage) has the form G^dl/dV ~ 
1 - B(T /Tk) 1 ^ g[AeV /k B T], where G is the T = 
linear-response conductance; A and B are constants that 
can be fitted to experiment. For x — > 0, g[x] — > const., so 
G^dl/dV ~ 1 - 5(T/T A -) 1/5 <?[0]; whereas g[x] ~ ex 1 ' 5 
(with c a constant) for x — > 00, yielding G^ 1 dl jdV ~ 
1 - cB(AeVlk B T K ) 1 ^ 25]. To lowest (quadratic) or- 
der in the tunneling matrix element between the impuri- 
ties and the lead [2fJ, the universal scaling function g[x] 
equals the exact function given in [11(b)]. Similar (and 111 
linear response, identical) behavior is expected in trans- 
port through triangular quantum-dot devices [23| . 

(c) Breaking of particle-hole symmetry. — This lowers 
the isospin SU(2) symmetry to the U(l) subgroup that 
conserves global charge 2I Z , while preserving the discrete 
S3 symmetry. The spectrum in Table IIIII is reclassi- 
fied by applying Eq. gj) to S[/«(2) 3 =) C/ W (l) x Z 3 . 
The most-relevant operators that become allowed in the 
low-energy Hamiltonian are marginal: the charge cur- 
rent operator 2I Z , which is exactly marginal and corre- 
sponds to a simple phase shift 16]; and a degenerate 
pair (s,I z ,Z 3 ,t z ,Z 8 ) = ((0)0,(0)0,(^)2/5,(0)0,(^)3/6) 
arising from Table IIIII row 5. The last two operators 
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are the boundary limits of the left- and right-moving 
bulk currents J L<R = 4>[ M [(T Z ) 2 - jl]ipL,R (Table |I| 
row 5) 113 ■ Jl,r generate a U(l) symmetry of the free- 
fcrmion bulk theory not preserved by the boundary con- 
dition. The boundary limits of such operators are exactly 
marginal |23j, consistent with NRG results in the pres- 
ence of particle-hole asymmetry Like a phase shift, 
the three exactly marginal deformations of the bound- 
ary conditions affect the FSS (and the boundary limits 
of Jl,r affect the T = zero-bias conductance), but not 
the operator spectrum in Table ITTT1 2^. Thus, the NFL 
fixed point and its signatures, including the ground-state 
degeneracy and power laws in the conductance, persist 
away from particle- hole symmetry (unlike, e.g., the NFL 
behavior of the two- impurity Kondo model |j|). 

(d) Other symmetry-breaking perturbations. — It can be 
deduced from Table ITTll that (i) spin-orbit coupling is rel- 
evant, with dimension 3/5, (ii) breaking of S3 symmetry 
(e.g., through distortion of the equilateral trian gula r im- 
purity geometry) is relevant with dimension 1/5 |28l | , (iii) 
spin-exchange anisotropy is irrelevant, (iv) a Zeeman field 
acting only on the impurity spins is exactly marginal, and 
(v) the coupling Jn in Eq. (J2J) is irrelevant. The impli- 
cations of these results will be discussed elsewhere [23J. 

In summary, we have found the exact low-energy be- 
havior of a non-Fermi-liquid phase arising from the in- 
terplay of magnetic frustration and Kondo physics in the 
three-impurity Kondo model. The phase is stable against 
particle-hole asymmetry, exchange anisotropy, and mag- 
netic fields. It should be detectable in tunneling into 
magnetic adatoms on metallic surfaces and in electrical 
transport through triangular quantum-dot devices. 
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